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We searched for scalar bottom quarks in 156 pb _1 of pp collisions at y/s = 1.96 TeV recorded 
by the CDF II experiment at the Tevatron. Scalar bottom quarks can be produced from gluino 
decays in R-parity conserving models of supersymmetry when the mass of the gluino exceeds that 
of the scalar bottom quark. Then, a scalar bottom quark can decay into a bottom quark and a 
neutralino. To search for this scenario, we investigated events with large missing transverse energy 
and at least three jets, two or more of which were identified as containing a secondary vertex from 
the hadronization of b quarks. We found four candidate events, where 2.6 ± 0.7 are expected from 
standard model processes, and placed 95% confidence level lower limits on gluino and scalar bottom 
quark masses of up to 280 and 240 GeV/c 2 , respectively. 

PACS numbers: 12.60.Jv 


Despite its extraordinary success, the standard model 
of elementary particles and interactions (SM) is in¬ 
complete since it does not explain the origin of elec- 


troweak symmetry breaking and the gauge hierarchy 
problem Jl]. One proposed extension of the SM, super- 
symmetry (SUSY) solves these problems. R-parity (2j 
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conserving SUSY models also provide a prime candi¬ 
date for the dark matter in the cosmos Q, the sta¬ 
ble lightest supersymmetric particle (LSP). SUSY is a 
space-time symmetry that relates particles of different 
spin by introducing a boson SUSY partner for each SM 
fermion and vice versa. For example, the left-handed 
and right-handed quarks have scalar partners denoted 
qx and qu which can mix to form scalar quarks (squarks) 
with mass eigenstates Several models jjf predict 

that this mixing can be substantial for the scalar bot¬ 
tom (sbottom), yielding a sbottom mass eigenstate ( 61 ) 
significantly lighter than other squarks. At the Tevatron 
center-of-mass energy of ^fs = 1.96 TeV, the predicted 
gluino ( g , the spin- 1/2 partner of the gluon) production 
cross section is almost an order of magnitude larger than 
that of a sbottom of the same mass |B|. Therefore, if 
sufficiently light, sbottom quarks could be copiously pro¬ 
duced through the decay of the gluino into a sbottom and 
a bottom quark since the gluino preferably decays into 
a squark-quark pair Q. A sbottom in the mass range 
accessible at the Tevatron is expected to decay predom¬ 
inantly into a bottom quark and a lightest neutralino 
which is often assumed to be the LSP. 

Previous searches for direct sbottom Q or gluino pro¬ 
duction [3] at the Tevatron placed upper limits on the 
masses of these particles. In this letter, we describe the 
first search for b\ from g decays in 156 pb _1 of data 
collected between July 2002 and September 2003 by the 
Collider Detector at Fermilab (CDF II) at the Tevatron. 
We assume R-parity conservation and, therefore, that g 
are produced in pairs. We consider a scenario where the 
branching fractions g — > bbi and 6 i —■» b\i are 100%. We 
assume a Xi mass of 60 GeV/c 2 which is above the cur¬ 
rent limits from LEP Q. The signature of our search is a 
final state with four 6 -jets from the hadronization of the 
6 -quarks, an imbalance in energy in the transverse plane 
to the beam (“missing transverse energy” or ifx) from 
the two LSPs escaping detection, and no isolated leptons 
with large transverse momentum. We expect that these 
events could be observed in inclusive three jet events with 
large After describing the experiment, we will intro¬ 
duce our event selection, major backgrounds, systematic 
uncertainties and results. 

CDF II is a multipurpose detector described in detail 
elsewhere E3- We use a cylindrical coordinate system 
around the beam axis in which 8 and (f> are the polar and 
azimuthal angles respectively. The pseudo-rapidity 77 is 
defined as 77 = — ln(tan(|)). The transverse momentum 
and energy are defined as px = P sin 8 and Ex = E sin 8 , 
respectively. The tracking system consists of a cylindrical 
open-cell drift chamber and silicon microstrip detectors 
in a 1.4 T magnetic field parallel to the beam axis. The 
silicon detectors provide precise tracking information for 
1771 <2 and are used to detect displaced vertices 0- 
The drift chamber m surrounds the silicon detectors 
and has maximum efficiency for |? 7 | < 1. The energies 


of electrons and jets are measured in calorimeters which 
cover 1 77 1 < 3.6. Muons are identified by drift chambers 
located outside the calorimeter volume which extend to 

M < 1-5. 

Jets are reconstructed from the energy depositions 
in the calorimeter cells using an iterative cone jet 
clustering algorithm [ll|, with a cone size of radius 
R = \J (A (j)) 2 + (A 77) 2 = 0.4. Energy corrections are ap¬ 
plied to account for effects that can distort the mea¬ 
sured jet energy such as non-linear calorimeter response, 
underlying events, or the position of the primary ver¬ 
tex. The ijx is defined by, iftx = — where 

n, is a unit vector perpendicular to the beam axis and 
pointing at the \ th calorimeter tower, and Elp the trans¬ 
verse energy therein. Candidate electrons must have a 
track associated with a cluster in the electromagnetic 
calorimeter with Ex > 10 GeV, an electromagnetic to 
hadronic energy ratio, and shower shape compatible with 
an electron 9j. Candidate muons are identified with 
px > 10 GeV/c tracks reconstructed in the tracking sys¬ 
tem which extrapolate to a track segment in the muon 
chambers, and leave an energy consistent with a mini¬ 
mum ionizing particle in the calorimeter :Qj. To increase 
the lepton efficiency (including r), candidate leptons are 
also identified as isolated tracks with px > 10 GeV/c. 

The data used for this search are selected by a three- 
level trigger system. The Level 1 trigger requires ijx > 
25 GeV , where the Tfx is estimated using trigger tow¬ 
ers with more than 1 GeV of total Ex- The Level 2 
trigger selects events with two jet clusters with Ex > 
10 GeV. At Level 3, the rfx is required to be larger 
than 35 GeV. After selections Q to remove accelerator- 
produced and detector-related background, and cosmic 
ray events, 418,304 events remain. The position of the 
primary vertex along the beam axis, Zo , is found by a 
weighted fit of each track’s z position at closest approach 
to the beam axis, and must satisfy |zo| < 60 cm. If there 
are multiple vertices in the event, we use the vertex with 
the largest J2(.Px) of associated tracks. The events are re¬ 
quired to have ifx > 35 GeV, and at least three jets with 
|? 7 | < 2 and Ex > 15 GeV. Events where the/^r is aligned 
with any of the three highest Ex jets are rejected by re¬ 
quiring Ad(^T, l-3jet) > where A(f> is the azimuthal 
angle between the ifix and each of the jets. The last se¬ 
lection criterion reduces backgrounds with tfx resulting 
from jet mismeasurements or from semileptonic 6 -decays. 
The long lifetime of B hadrons yields secondary vertices 
which have a large positive decay distance L xy [5] rela¬ 
tive to the primary vertex. We require the events to have 
at least one jet identified (inclusive single 6 -tagged) as a 
6 -jet by the CDF 6 -tagging algorithm [ljj. After apply¬ 
ing these selections, 332 events are observed in the data. 

Dominant SM backgrounds are multijet, top, and 
IU/Z+jets/diboson (electroweak) production. The mul- 
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FIG. 1: Distributions of A 4>(I^t, 1 s ^ jet), the azimuthal angle between the I/t and the highest Et jet in the event, if it is 
6 -tagged in control region CO (left), Cl (middle), and C2 (right). 


TABLE I: Comparison of the total number of expected and observed inclusive single 6-tagged events in the control regions. 
Only statistical uncertainties are quoted. 


Control region 

CO 

Cl 

C2 


35 - 50 GeV 

35 - 50 GeV 

> 50 GeV 

isolated high-pr lepton 

required 

vetoed 

required 

W /Z+jets/diboson 

3.9 ±0.8 

10.9 ± 1.2 

9.6 ± 1.2 

Top 

11.7 ± 0.2 

8.2 ±0.1 

35.0 ±0.3 

Multijet 

19.1 ± 4.1 

128.9 ± 17.3 

10.9 ±4.5 

Total predicted 

34.7 ±4.2 

148.0 ± 17.3 

55.5 ±4.7 

Observed 

36 

121 

63 


tijet background is composed of a heavy flavor (HF) com¬ 
ponent and light jets misidentified by the 6 -tagging algo¬ 
rithm as HF jets. The latter are estimated from data us¬ 
ing secondary vertices with negative L xy . The PYTHIA 
event generator 0 is used to estimate the HF and top 
background. The top production is normalized to theo¬ 
retical predictions [17j while the HF contribution is scaled 
to the data 0. The electroweak background is gener¬ 
ated in ALPGEN 0 , with the cross section at next-to- 
leading order (NLO) calculated by the MCFM 0 pro¬ 
gram, and the partons are showered in the HERWIG 0 
program. The CTEQ5L 0] parton distribution func¬ 
tions are used in the generation of all background sam¬ 
ples. Events are passed through the standard simula¬ 
tion j|22j of the CDF II detector and weighted by the 
probability that they would pass the trigger as deter¬ 
mined in independent data samples. 

To avoid potential bias when searching for new physics, 
we tested the SM background in control regions that are 
defined a priori. The data were divided into four regions 
having different ifx and high-pp isolated lepton multiplic¬ 
ity. The signal region contained events with large l/p and 
no high-pr isolated leptons. The three control regions 
used to check the SM prediction, denoted as CO, Cl, and 
C2, are defined in Tabled CO and Cl are sensitive to the 
predictions of the multijet background. Control region 
C2 with high-pT isolated leptons and I/t > 50 GeV is 
used to check the top and electroweak backgrounds. Pre¬ 


dicted total number of events and distributions of kine¬ 
matic variables such as the jet Et , the track multiplicity 
and the 1 /t have been studied and found to be in agree¬ 
ment with observations Q in the control regions. As an 
example, Fig.UJshows A0(^r, 1 s ^ jet), the azimuthal an¬ 
gle between the ifir and the highest Et jet in the event 
(l st jet) if this jet is 6 -tagged. The distributions show 
good agreement between data and prediction. Table d 
summarizes the background contributions to the num¬ 
ber of expected inclusive single 6 -tagged events and the 
observed events in the control regions. 

We optimize the sensitivity to sbottom produc¬ 
tion from gluino decays by maximizing the ratio of 
Nsusy/VNsm where Nsusy and N$m are the sbot¬ 
tom and SM background events expected in the sig¬ 
nal region [§.], respectively. The signal predictions 
are obtained by computing the acceptance using the 
ISAJET |0 event generator normalized to the NLO pro¬ 
duction cross section determined with PROSPINO f| 
and the CTEQ5M |0 parton distribution functions. 
Figure O shows the I^t spectrum for events with a single 
6 -tagged jet (exclusive single 6 -tagged) and at least two fa- 
tagged jets in the event (inclusive double 6 -tagged) when 
high-pT isolated leptons are vetoed. The best sensitivity 
is achieved by requiring inclusive double 6 -tagged events 
with IpT > 80 GeV. The double 6 -tag requirement is 
very effective at suppressing the background. Exclusive 
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single tagged events were used to cross-check the results. 
Since the efficiency of tagging a 6-jet is about 30% 0, 
and there are four 6-jets in the signal final state, the 
signal acceptance for single and double tag events are 
similar and vary between 5% and 12% as a function of 
the sbottom and the gluino masses. 

The systematic uncertainty on the signal and the back¬ 
ground predictions is studied in detail {$]. Correlated 
and uncorrelated uncertainties are evaluated separately 
and combined in quadrature. Dominant uncertainties 
that affect both the background prediction and the sig¬ 
nal are the jet energy scale (25% for multijet and elec- 
troweak, 14.5% for top, 10% for the signal), 6-tagging 
efficiency (12%), and the luminosity (6%). Uncorrelated 
systematic uncertainties on the background predictions 
are dominated by the heavy flavor MC scale factor (20%), 
the misidentified 6-tag rate (16% for light flavor multi- 
jets), the top cross section (15%), and the electroweak 
cross section (11.5%). The total systematic uncertainties 
are 35%, 31%, and 25% for the multijet, electroweak and 
top background predictions, respectively and 18% on the 
signal. 

The signal region is analyzed after all the background 
predictions and selection cuts are determined. Requiring 
at least two 6-tagged jets we observed four events, where 
2.6 ± 0.7 are expected, as summarized in Table HU We 
observe 21 exclusive single 6-tagged events, which is in 
agreement with SM background expectations of 16.3±3.6 
events. 

TABLE II: Number of expected and observed events in the 
signal region. Correlated and uncorrelated uncertainties in 
the total predicted background were treated separately and 
combined as described in Q. 


Process 

Excl. Single 6-Tag 

Inch Double 6-Tag 

W /Z+jets/diboson 

5.6 ±1.8 

0.6 ±0.3 

Top 

6.1 ±1.5 

1.9 ±0.5 

Multijet 

4.6 ± 1.7 

0.2 ±0.1 

Total predicted 

16.3 ±3.6 

2.6 ±0.7 

Observed 

21 

4 


Since no evidence for gluino pair production with 
g —> 66i is found, a cross section upper limit at 95% confi¬ 
dence level (C.L.) is computed and an exclusion limit set 
using the Bayesian likelihood method with flat prior 
probability for the signal cross section and Gaussian pri¬ 
ors for the uncertainties on acceptance and backgrounds. 
An example cross section upper limit is shown in Fig. |31 
The contour exclusion limit shown in Fig. |3j is computed 
using the inclusive double 6-tag analysis. The mass 
has a minor impact (9| on this limit, as long as the Xi is 
sufficiently lighter than the 6i, so that its decay produces 
a taggable 6-jet. 

In conclusion, we have searched for sbottom quarks 
from gluino decays in 156 pb -1 of CDF Run II data. We 



P T [° eV ] 



P T [Ge y ] 

FIG. 2: I^t spectrum after vetoing events with high-py iso¬ 
lated leptons, for single 6-tagged events (top) and double 6- 
tagged events (bottom). 


observe four inclusive double 6-tagged candidate events, 
which is in agreement with SM background expectations 
of 2.6 ±0.7 events. No evidence for sbottom quarks from 
gluino decays is observed and we exclude a significant re¬ 
gion in the gluino and sbottom mass plane at 95% confi¬ 
dence level. Gluino masses below 270 GeV/c 2 for Am = 
m(g) — m(6i) > 6 GeV/c 2 and m(6i) < 220 GeV/c 2 are 
excluded. 
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FIG. 3: The 95% C.L. cross section upper limit as a function 
of the gluino mass for a 60 GeV/c 2 mass difference between 
the gluino and sbottom (top). The 95% C.L. limit curve in the 
m(g)-m(bi) plane for inclusive double 6 -tagged events (bot¬ 
tom). 
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